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Abstract 

We try to explain the spatial variation of primordial deuterium suggested 
by some observations by varying leptonic chemical potentials. The variation of 
the latter may take place in some scenarios of leptogenesis. The model predicts 
a large mass fraction of 4 He (35-60%) and 1 Li (up to 10~ 9 ) in deuterium- 
rich regions. Because of lepton family symmetry, the angular variations of 
cosmic microwave background radiation can be sufficiently small although still 
observable in future measurements. 

Recently several groups [1-6] have reported measurements of the deuterium abun- 
dance in Lyman-limit absorption line systems with red-shifts 0.48 < z < 3.5 on the 
line of sight to quasars; these are believed to give essentially the primordial value. 
Surprisingly some groups have claimed a high value, D/H ~ 2 ■ 10~ 4 on the basis 
of ground-based data taken with the Keck telescope, but this result is now thought 
to be due to various errors and the best value available from two "clean" sys- 
tems is 3 ■ 10~ 5 H. However, Webb et al HI report a high deuterium abundance, 
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D/H ~ 2 • 10~ 4 , in an apparently clean system with z = 0.7 observed with the Hub- 
ble Space Telescope, as well as a low one in another system with z = 0.5, raising the 
possibility that there might be real spatial variations in primordial D/H. 

If the effect is indeed real (which it is perhaps too early to judge), its significance 
is difficult to overestimate. It would strongly change our approach to primordial nu- 
cleosynthesis and possibly to the physics of the early universe. A possible variation 
of the light element abundances was in fact considered in ref. (see also ||), where 
a model of leptogenesis was considered which, first, gave a large lepton asymmetry, 
which could even be close to or larger than 1, and, second, this asymmetry might 
strongly change on astronomically large scales, It,. The magnitude of the latter de- 
pends on the unknown parameters of the model and can easily be in the mega-giga 
parsec range. The model is based on the Affleck-Dine || scenario of baryogenesis but 
in contrast to the original one it gives rise to a large (and varying) lepton asymmetry 
and to a small baryonic one. Recently a similar model of generation of large (but 



not varying) lepton asymmetry was considered in ref. [I0|. In what follows we will 
not discuss the details of the model but confine ourselves to a more phenomenological 
level, namely we simply assume that there exists a mechanism which created large 
leptonic asymmetries of order unity (electronic, muonic and/or tauonic) which vary 
by 100% over the distance 1^. A possible early universe scenario which would give 
rise to such varying and large leptonic asymmetries will be considered elsewhere. At 
the moment we put a less ambitious question: whether it is possible to describe the 
suspected spatial variation of primordial deuterium by varying chemical potentials of 
neutrinos without conflict with the existing astronomical data and what predictions 
can be made in such a model which can be tested in future observations. 

To explain the claimed variation of 2 if at z — 0.7 the magnitude of II in terms of 
present-day units must be smaller than or close to one gigaparsec. The lower bound on 
this scale II may be much smaller. It can in principle be determined by measurements 
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of the abundances of light elements at large distances in our neighborhood, say, z > 
0.05. It would be interesting if the scale II coincides with the 140/h Mpc scale 
observed in the large scale structure of the universe [IT], IT2 . 



Another simple possibility to explain a varying abundance of deuterium is to 
assume that the baryon-to-photon ratio in the universe varies as a function of position. 
This idea was studied in ref. |13| where it was shown that the necessary large scale 



isocurvature perturbations are excluded by the smallness of angular fluctuations of 
the cosmic microwave background radiation (CMB). A similar criticism is applicable 
at first sight to the model with varying lepton asymmetry. Indeed, it can be easily 
checked that the necessary change in chemical potential of electron neutrinos £„ e 
should be close to —1 to explain the possibly observed variation of deuterium by 
roughly an order of magnitude. Such a change in £ Ve would induce a variation in 
total energy density during the RD stage at a per cent level, which is excluded by the 
smoothness of CMB. However, this objection can be avoided if there is a conspiracy 
between different leptonic chemical potentials such that in different spatial regions 
they have the same values but with interchange of electronic, muonic and/or tauonic 
chemical potentials. In other words we assume that in a particular spatial region the 
three neutrino chemical potentials have the values 

= 7]- (1) 

Then in another spatial region they should have the same values but with an arbitrary 
interchange of e, //, and r. This would ensure the same energy density at different 
space points and small angular variations of CMB. In fact the perturbations in CMB 
induced in this way would be non- vanishing and close to existing observations. We will 
discuss them below. Since the abundances of light elements are much more sensitive 
to the magnitude of the electron neutrino chemical potential than to those of muon 
and tauon neutrinos, the variation of ^ Ve (accompanied by corresponding variations 



of and ^ Vt ) would lead to a strong variation in the abundance of deuterium and 
other light elements. 

The equality of, say, £„ e at one space point to at another point looks like 
a very unnatural fine-tuning but this is not so. The present theory of elementary 
particles is believed to be symmetric with respect to interchange of three families of 
leptons. In the Affleck-Dine type scenario of generation of charge asymmetry, the 
latter is generated owing to the formation of baryonic (as in the original version) 
or leptonic charge condensates along the so called flat direction in the potential of 
a scalar field which possesses corresponding charges. It is rather natural to assume 
that the potential respects the symmetry between different lepton families. So if a 
flat direction corresponds to a scalar field with the combination of leptonic charges 
[a,/3, 7], then there must be flat directions with the same values of the leptonic 
charges but interchanged with respect to e, fi, and r. In such a model there would 
be regions with different values of leptonic chemical potentials which are obtained by 
transmutations of the original ones in (JlJ) . 

The symmetry between lepton families is broken at low energy by the masses of 
charged leptons. So one may expect that there could be significant fluctuations of the 
cosmic energy density when the temperature is close to the mass of the charged tau- 
lepton (m = 1777 MeV) or to that of the muon (m = 106 MeV). We will see below that 
this is not the case. Another potential danger for a model of this kind is the variation 
of the energy density associated with the energy of the potential wall between the 
valleys (flat directions) with different leptonic charges. These domain walls however 
disappear because all scalar fields (pi possessing different leptonic charges evolve down 
to the same origin of the potential where they all vanish, (pi = 0. Some remnants of 
the energy distortion remain but they have an energy density much smaller than that 
of the original domain walls. 

Let us turn now to the astro-phenomenology of the model concerning the change 
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in the light element abundances. It is straightforward and simple to play with the 



standard nucleosynthesis code |L4| to study the influence of different leptonic chemical 
potentials on the output of light elements and in Table 1 we present a number of sam- 
ple calculations from which we can draw some combinations that would give a small 
amount of 2 H in our neighborhood (and a larger one in deuterium- rich regions). Con- 
siderations of Galactic chemical evolution [H| permit us to infer that the abundance 
of primordial deuterium in nearby regions where A He is also measured is close to the 
low values determined at high red-shift; we take as the best estimates for both these 
regions and our neighborhood D/H = (3.2 ± 0.8) • 1(T 5 and R( 4 He) = 0.24 ± 0.01, 
which are well fitted in the case of no neutrino degeneracy for a baryon/photon ratio 
7710 = 5 ± 1. An adequate fit is also obtained for the same 77 if we take the combi- 
nation [£i/ e ,£i/ ] = [0, — 1], which gives for the "mirror" region with [^ e ,^ ] = [—1,0] 
a substantially higher deuterium abundance D/H = 8.5 • 10~ 5 . This combination is 
not necessarily the best possible fit to the data, be the but it seems too early to look 
for this, bearing in mind that the observational data may change. It is worth noting 
that if two (or all three) £'s are permitted to vary, the nucleosynthesis limits (for a 
recent reference see e.g. |[L6| ) would be invalidated. 

One can see from Table 1 that the data somewhat resist the proposed explanation. 
It would help if there is more deuterium in our neighborhood, ~ 5 • 10~ 5 , and/or less 
in the deuterium-rich regions, ~ 10~ 4 . We did not try to use large values of chemical 
potentials because of possible problems with smoothness of CMB temperature (see 
below). The agreement with observations can be made better if all three chemical 
potentials could be adjusted as free parameters. In Table 2 we present the abundances 
of light elements for the ad hoc choice [£„ e , £„ , £ Vt ] = [—1, 0.1, 1] for r] W = 4 and 5. The 
last line may not be reliable because the program fails to converge. It is noteworthy 
that it is possible to have, besides high 2 H and 4 He regions, regions with normal 
deuterium and low helium-4. 
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To describe simultaneously the suspected deuterium content in the rich regions 
and that in our neighborhood we need 77x0 = 5 — 6 and £ Ue =0 — 0.1 in the poor regions 
and £ Ue ~ — 1.4 in the rich regions. In this case it is possible to get D/H as large as 
17- 10~ 5 in rich regions. The necessary value of £ is rather high and it would be easier 
for the model if the deuterium fraction in rich regions would be around 10 • 10~ 5 . 

A generic feature of our model is that simultaneously with high deuterium a high 
mass fraction of helium-4 is predicted. It is at least 30-35% or may be even above 50%. 
It is an interesting question what is the observational upper bound on the abundance 
of A He far away from us. All direct measurements of A He known to us were done at 
most at z = 0.045 corresponding to a distance of 140/i _1 Mpc||17||. A very large mass 
fraction of 4 ife can possibly be excluded with the help of star and galaxy evolution. 
Stars should be brighter and have a shorter life-time. All data indicate that distant 
objects (including quasars) have more or less normal chemical content. Still we do 
not know what is the permitted mass fraction of helium-4 which does not contradict 
the data. This would be the subject of a separate study. Presumably 35% of 4 He in 
some distant parts of the universe is not excluded. As for much higher values, we do 
not have an answer now. 

A very sensitive indicator of any inhomogeneities in the universe is the cosmic 
microwave background. The bearers of electronic, muonic, and tauonic chemical 
potentials have different masses: though different neutrinos are most probably very 
light or even massless, so that their contribution to the energy density is the same, the 
masses of the corresponding charged leptons are very much different and this could 
be potentially dangerous for the model. This is not the case, however, as can be seen 
from the following considerations. Let us assume for simplicity that there are only two 
lepton families, electronic and muonic. Let us assume also that the primeval plasma 
has nonzero electronic and muonic charge densities, D e and D^. Of course the plasma 
is electrically neutral. Thermal equilibrium in the plasma is fulfilled with a very good 
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accuracy (at least for temperatures above 3 MeV, when neutrinos decouple). The 
distributions of different particles are given by the normal Fermi (or Bose) functions 
with nonzero chemical potentials which permit to have nonzero D e and D^. Due to 
reactions e~ + v e <-> \l~ + and similar (crossed) ones, the following relation between 
chemical potentials must be fulfilled in thermal equilibrium: 

£e £,v e = C/ti bp (2) 

There is also the condition of electric neutrality of the plasma: 

5n e + 5n^ = (3) 
and the expressions for electronic and muonic charge densities: 

Sn e + 5n Ve = D e (4) 

and 

5n^ + 5n^ = (5) 

where 5n a = n a — n a is the difference in number densities of particles and antiparticles 
with 

Ua = / ~P~i TWl^ — FT ( 6 ) 

J 1 + exp(E/T - ia) 

The number density of antiparticles is given by the same expression with the opposite 
sign of £ . 

One can see from the symmetry property of the system of equations (§-[|) that 
for any solution corresponding to the set \D e ,D^\ = [a, (3\ there exists the mirror 
solution corresponding to the set [D e , DJ\ = [j3, a] which can be constructed from 
the original solution by the substitution: 5n e «-» —5n e , Sn^ <-> — Sn^ (correspondingly 
£e,n —£,e,n) and 8n Ve <-» 5n Vll . Evidently the energy densities of both solutions are 
the same. 
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There are some other possible inhomogeneities in the energy density that could 
be either dangerous for the model or observable in CMB. The first and most evident 
one is related to the binding energy of 4 He, which is 7 MeV per nucleon. Since the 
mass fraction of 4 ife may change by a factor of 2 in deuterium- (and helium-) rich 
regions (from 25% to more than 50%), this means that the variation in baryonic 
energy density may be as large as 2 • 10~ 3 . The contribution of baryons to the total 
energy density is given by Qb = 3%(?7/4)(0.65//i) 2 , so the relative density fluctuations 
are at most 

^ = 6-l(T 5 (77/4)(0.65//i) 2 (7) 

Ptot 



To estimate the fluctuations in CMB temperature we can use the results of ref. |fL3 
where similar isocurvature density perturbations, but with amplitude (2 ■ lO -3 ) -1 = 
500 x larger, were considered. According to their results normalized to our smaller 
perturbations 

- = l0- 5 (^) 2 (8) 

where Ao = c/H = 3Gpc//i. So the anisotropy induced by the Sachs- Wolfe effect 
would be below the observational bounds for scales above ~ SOOh^ 1 Mpc. If D/H 
in the rich regions is about 10~ 4 (and not 2 ■ 10~ 4 ), then the variation of helium-4 
could be smaller. Correspondingly smaller density perturbations would be induced. 
In this case smaller scales in fluctuations of CMB temperature, down to ~ 150/t -1 
Mpc, would be permitted. However, on scales below 2° (or below 200 Mpc), the 
result (^|) would not be valid. These scales are dominated by Doppler shift across the 
fluctuations at the surface of last scattering [PJ. The measurements on small scales 
permit possibly 5T/T = 3 ■ 10~ 5 , which could also be compatible with this model. 
Such fluctuations may be observed in the future MAP or PLANCK missions or with 
balloons. 

There is another effect which is more subtle theoretically but which could also 
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give rise to similar fluctuations in 5T/T. The energy densities of electron and muon 



neutrinos are known [[L9| , |20|| to be different owing to the following effect. After 
neutrinos decoupled from the primeval plasma, which roughly took place at T = 2 
MeV for electron neutrinos and at T = 3 MeV for muonic and tauonic ones, the 
temperatures of electrons and photons became somewhat different from the neutrino 
temperature owing to heating of the electromagnetic component of the plasma by 
e + e~-annihilation into photons. Because of this temperature difference and due to 
residual e + e~-annihilation into vv, the usually assumed equilibrium neutrino distribu- 
tions became slightly distorted. The nonequilibrium correction to the energy density 
of electron neutrinos in the standard model is approximately ||21||: 



ApJ Pu « 0.9% (9) 

and the distortion of the energy density of muon and tauon neutrinos is 

Ap u Jp u = A Pu J Pu « 0.4% (10) 

(closely similar results are obtained in ref. [0). The difference between v e and v^ T 
is related to a greater efficiency of the process e + e~ — > vv due to the presence of 
charged current interactions only for v e . Now because of nonzero leptonic chemical 
potentials these results would slightly change. They would remain the same in the 
Boltzmann approximation because the probability of e + e~-annihilation into vv does 
not depend on the chemical potential of neutrinos in the case of Boltzmann statistics. 
Typically corrections due to Fermi statistics are about 10%. So the relative efficiency 
of annihilation due to nonzero chemical potentials becomes smaller by approximately 
0.1 [cosh (£) — 1]. This expression is true for relatively small £, £ < 1; for larger £ it is 
changed to a power law. 

To get an estimate of the magnitude of the density fluctuations due to variation 
of chemical potentials let us assume that in our neighborhood the chemical potentials 
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have the values £ Ve = £ Ut = and £„ = — 1 and in the deuterium- rich region they are 
£uu — iu T = and £„ e = —1. Thus the relative energy density of neutrinos changes by 

Ml = S (^ + + = (0.9%-0.4%)-0.1(coshf-l) ~ 2.5-10' 4 (11) 

Pu \ Pv Pv Pu J 

Keeping in mind that one neutrino species contributes 10-20 % to the total energy 
density during the RD stage, we find that the relative density fluctuations of neutrinos 
due to variation of chemical potentials are approximately bp v j ptot ~ 5 ■ 10~ 5 . In fact 
the fluctuations of the total energy density are very much smaller than that because 
the increase in p v is accompanied by a similar decrease in the energy density of 
photons and e^. Thus the phenomenon we discuss gives rise to a rather peculiar 
perturbation: the variation of the total energy density is negligibly small but the 
radiation temperature varies between different spatial points. 



As was calculated in ref. [21| the photon temperature drops in comparison with 
the standard one by 10~ 3 , due to the above mentioned transfer of energy from the 
electromagnetic component of the plasma to neutrinos. This change of temperature 
should be proportional in a crude approximation to the above mentioned change of 
neutrino energy: 

^O.lMl. 1(r 3 (12) 

T p v 

Now if chemical potentials are not spatially constant, the quantity Ap^" 1 ' / p v would 
vary together with the chemical potentials at different points. Its variation is given 
by eq. (]TT|) . Accordingly the variation of the photon temperature due to this effect is 

— = 5 1 0.1-^— w 5 • 10- 5 (cosh£ - 1) ~ 2.5 • 10" 5 (13) 
T V Pv J 

which is close to the observational bounds. 

These are of course very crude estimates. The real result should be somewhat 

smaller. An account of inverse annihilation, z/ e z/ e — > e~e + and of elastic z/e-scattering 

results in a smoothing down of the spectral distortion. An estimate of these "inverse" 
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effects, made along the lines of semi-analytical estimates of ref. [19], diminishes the 
temperature change by a factor of roughly 2/3. 

The magnitude of temperature fluctuations depends in particular on the unknown 
values of the £'s. For example for |£| = 0.7 the effect would be twice smaller than 
for |£| = 1, while for |£| = 1.4 it is twice bigger. We take |£| = 1.4 as an upper 
limit for the magnitude of possible variations of chemical potentials. To be on the 
safe side we possibly need somewhat smaller £'s and correspondingly a fraction of 
deuterium in the rich regions of about 10~ 4 . A rigorous calculation of the effect is a 
straightforward but formidable numerical problem. It seems premature to do that at 
this stage. However, we will have to do the calculations if the effect of spatial variation 
in deuterium abundances is confirmed and the predicted variation of helium-4 is either 
found or not ruled out. 

To conclude, we try to explain a spatial variation of primordial deuterium, that has 
perhaps been observed, by varying leptonic chemical potentials. The model could be 
confirmed (or rejected) by looking for a very large mass fraction of primordial helium- 
4 in deuterium-rich regions but this is not a practical possibility in the context of data 
available now or in the near future. A more promising test seems to be possible from 
the theory of stellar evolution with a high mass fraction of 4 He. The hypothesis also 
predicts larger abundances of other light elements in these regions, e.g. 7 Li should 
be at the level of 10 -9 . There might be also regions with normal deuterium and low 
helium-4. If there is a family symmetry which ensures permutational symmetry for 
different chemical potentials, a very large distortion of CMB isotropy can be avoided, 
but there still remain nonzero 5T/T fluctuations which can be detected if the higher 
abundances of deuterium and other elements (in particular 4 He) exist. It is a curious 
coincidence that the theory of large scale structure formation may possibly favor 
neutrino chemical potentials |23|] close to those that are needed in our model. 
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Table 1: Abundances of light elements for different values of the baryon number 
density, r\ = 10 10 71b/ti 7 and neutrino chemical potentials £„ o . 
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Table 2: Abundances of light elements for rj = 10 10 n£/n 7 = 4, 5 and nonzero values 
of all three neutrino chemical potentials £ Va . 
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